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^v^ i We present an algoritlim for calculating tlie complete data on an event 

K*" I liorizon which constitute the necessary input for characteristic evolution of 

Qv^ ' the exterior spacetime. We apply this algorithm to study the intrinsic 

^D ' and extrinsic geometry of a binary black hole event horizon, constructing 
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a sequence of binary black hole event horizons which approaches a single 



^^ , Schwarzschild black hole horizon as a limiting case. The linear perturbation 

^-^ ' of the Schwarzschild horizon provides global insight into the close limit for 



binary black holes, in which the individual holes have joined in the infinite 
past. In general there is a division of the horizon into interior and exterior 



bJO' regions, analogous to the division of the Schwarzschild horizon by the r = 2M 

bifurcation sphere. In passing from the perturbative to the strongly nonlinear 
regime there is a transition in which the individual black holes persist in the 



$H ' exterior portion of the horizon. The algorithm is intended to provide the data 

sets for production of a catalog of nonlinear post-merger wave forms using the 
PITT null code. 
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I. INTRODUCTION 

In previous work, we developed a model which generates the intrinsic null geometry of 
an event horizon with the "pair-of-pants" structure characteristic of a binary black hole 
merger |jl|,0]. In this paper, we extend this approach to determine all extrinsic curvature 
properties of such horizons, thus providing a complete stand-alone description of the event 
horizon of binary black holes. We apply this work to study the event horizon of a head-on 
collision of black holes using a sequence of models which embraces not only the perturbative 
regime of the close approximation 0, where the merger takes place in the distant past, 
but also includes the highly nonlinear regime. In the perturbative regime the individual 
black holes merge in an interior region of the horizon, corresponding to the region of the 
Schwarzschild event horizon lying inside the r = 2M bifurcation sphere. But we show how 
dramatically nonlinear effects can push the merger into the exterior portion of the horizon. 

Beyond the investigation of the horizon geometry of a black hole collision, the major 
motivation for this work is to provide the null data necessary to compute the emitted gravi- 
tational wave by means of a characteristic evolution of the exterior spacetime. In the Cauchy 
problem, the necessary data on a spacelike hypersurface are the intrinsic metric and extrinsic 
curvature, subject to constraints. On a null hypersurface, such as an event horizon, the situ- 
ation is quite different. The necessary null data consist of the conformal part of the intrinsic 
(degenerate) metric, which can be given freely as a function of the affine parameter. Then 
the surface area of the horizon is determined, via an ordinary differential equation along the 
null rays (the Raychaudhuri equation), in terms of an integration constant supplied by the 
mass of the final black hole. Similarly, all extrinsic curvature components of the horizon are 
determined by ordinary differential equations in terms of integration constants supplied by 
the final black hole. Whereas in principle the intrinsic conformal geometry is the only data 
on the horizon necessary for the characteristic initial value problem, in practice the surface 
area and extrinsic curvature are essential to supply the start-up data for the implementation 
of a characteristic evolution code, such as the PITT null code 0,^. The main results of 
this paper concern the understanding of the nonlinear nature of the underlying ordinary 
differential equations from geometrical, physical and numerical points of view. 

Given the intrinsic geometry and extrinsic curvature of the horizon, the strategy be- 
hind the characteristic approach to the computation of the emitted wave has been outlined 
elsewhere 0. It consists of two evolution stages based upon the double null initial value 
problem. Referring to Fig. 0, the necessity of two stages results from the disconnected 
nature of the two null hypersurfaces on which boundary conditions must be satisfied: (i) 
the event horizon 7i+, where binary black hole data is prescribed, and (ii) past null infinity 
X~, where ingoing radiation must be absent, at least in the late stage to the future of the 
hypersurface F which is decoupled from the formation of the individual black holes. Rather 
than directly attempting to solve this mixed version of a characteristic initial value problem, 
the stage I evolution is based upon two intersecting null hypersurfaces consisting of the event 
horizon, where binary black hole data is prescribed, and an ingoing null hypersurface J^ 
approximating future null infinity X"*", where data corresponding to no outgoing radiation 
is prescribed. The characteristic evolution then proceeds backward in time along ingoing 
null hypersurfaces extending to I~ to determine the spacetime exterior to an ingoing null 
hypersurface F. As indicated in the figure, the evolution terminates at F because the horizon 



splits into two individual black holes. This stage I solution is the advanced solution to the 
problem in the sense that radiation from 1~ is absorbed by the black holes but no outgoing 
wave is radiated to X+. Stage II provides the retarded solution, where outgoing but no ingo- 
ing radiation is present outside F, by running forward in time a double null evolution based 
upon the intersecting null hypersurfaces F, where the stage I data is prescribed, and an 
outgoing null hypersurface J~ approximating X~, where data corresponding to no ingoing 
radiation is prescribed to the future of F. The stage II evolution produces the retarded 
solution for the spacetime outside the world tube F. The approach is a nonlinear version of 
the standard method of determining the retarded solution ^ret of the linear wave equation 
□ $ = 5*, with source S, by first finding the advanced solution ^adv and then superimposing 
the source-free solution ^ret — ^adv- In the nonlinear case, where standard Green function 
techniques cannot be used to define retarded and advanced solutions, they can be defined by 
requiring the absence of radiation at X"*" (retarded) or X~ (advanced). In the present case, 
the data on the world tube F represents the interior source which has led to the formation 
of the individual black holes. (In a normal physical context, the source consists of two stars 
undergoing gravitational collapse but in a purely vacuum scenario imploding gravitational 
waves can play the role of the matter.) The justification of this two stage approach is that it 
reduces to the standard linear method in the close approximation where the geometry can be 
regarded as a perturbation of a Schwarzschild background. A purely perturbative character- 
istic treatment of the close approximation using the characteristic approach has been carried 
out in a separate study and the advanced solution has been successfully computed 0. 

We intend to carry out the "backward in time" evolution using the PITT null code [^|1 , 
which is based upon the Bondi-Sachs version of the characteristic initial value problem [§|,3. 
The present code is designed to evolve forward in time along a foliation of spacetime by 
outgoing null hypersurfaces. In this paper, in order to apply it to the backward in time 
evolution of a black hole horizon, we describe the merger of two black holes in the time 
reversed scenario of a white hole horizon, in which the black hole merger is now represented 
by the fission of a white hole. The post-merger era of the black hole horizon then corresponds 
to the pre-fission era of the white hole; and the proposed backward in time evolution of the 
black hole horizon to determine the exterior spacetime corresponds to a forward in time 
evolution of the white hole horizon. Our algorithm provides the complete white hole data 
necessary to carry out this evolution. 

The specific application in this paper is to the axisymmetric head-on collision of two 
equal mass black holes. However, the algorithm is capable of generating intrinsic geome- 
try and extrinsic curvature of an arbitrary event horizon, including the case of inspiraling 
binary black holes of non-equal mass. In previous studies of the intrinsic geometry of non- 
axisymmetric horizons, the approach has revealed new features of the generic collision of two 
black holes, such as an intermediate toroidal phase which precedes the merger [^ . Here we 
apply the algorithm to gain new insight into the global behavior of the extrinsic curvature 
properties of the head-on collision. 

The material in the paper divides into two types: (i) general features of a binary horizon 
and (ii) the calculation of the null data on the horizon necessary to initiate the PITT code. 
Material of the first type, which describes how a binary horizon deviates from the close 
approximation in the nonlinear regime, appears in Sees. [II A| and |T| - ^ and can be read 



independently. This material is described most conveniently in terms of Sachs coordinates, 



introduced in Sec. [II A|. Material of the second type involves the transformation from Sachs 



to Bondi-Sachs coordinates, introduced in Sec. II B| . This material supplies all data necessary 



to compute, using the PITT code, the fully nonlinear merger to ringdown wave form, which 
will be the subject of future work. The details of implementing the horizon data algorithm 
as a finite difference code are given in the Appendixes. An independent code designed for 
an axisymmetric horizon is also described and has been used as an independent test of the 
full code. 

We retain the conventions of our previous papers [||,^^3, with only minor changes in 
notation where noted in the text. For brevity, we frequently use the notation f^ = d^f to 
denote partial derivatives and / = duf to denote retarded time derivatives. 

II. CHARACTERISTIC DATA ON A HORIZON 

A. Horizon structure 

The evolution of the exterior spacetime by the PITT code proceeds along a family of 
outgoing null hypersurfaces. The characteristic initial value problem for the evolution re- 
quires an inner boundary condition which can be set either on a timelike world tube or, as a 
limiting case, on a null world tube. Here we choose the inner boundary to be the null world 
tube representing a white hole horizon Ti. The white hole horizon pinches off in the future 
where its generators either caustic or cross each other (such as at the vertex of a null cone). 
As illustrated in Fig. ^, we introduce (i) an afiine null coordinate u along the generators of 
7Y, which foliates 7i into cross sections Su and labels the corresponding outgoing null hy- 
persurfaces Ju emanating from the foliation; (ii) angular coordinates x^ which are constant 
both along the generators of Ti and along the outgoing rays and (iii) an afiine parameter A 
along the outgoing rays normalized by V"mVq,A = —1, with A = on ?i. In the resulting 
x" = (m, X,x^) coordinates, the metric takes the form 

ds"^ = -{W - gABW^W^)du^ - 2dud\ - 2gABW^dudx'^ + gAsdx^dx^ . (2.1) 

The contravariant components are given by g^^ = —1, g^^ = —W"^, g^^ = W and g^^gsc = 
6q. In addition, we set g^B = r'^hAs, where det^hAs) = det{qAB) = <l{x^), where qab is 
some standard choice of unit sphere metric. 

These coordinates were first introduced by Sachs to formulate the double-null character- 



istic initial value problem |T^. The Bondi-Sachs coordinate system p,^ used in the PITT 
code differs by the use of a surface area coordinate r along the outgoing cones rather than 
the afiine parameter A. However, because the horizon is not a surface of constant r, except 



in the special case of an "isolated horizon" [|Tl|, it is advantageous to first determine the 



necessary data in terms of an afiine parameter and later transform to the r-coordinate. 

The requirement that Ti be null implies that PF = on ?i. There is gauge freedom on 
7i that we fix by choosing the shift so that du is tangent to the generators, implying that 
W^ = on 7i; and by choosing the lapse so that u is an afiine parameter, implying that 
d\W = on 7"^. We adopt these choices throughout the paper, and our results generally 
hold only on 7f (A = 0) and when these conditions are satisfied. Later, in Sec. |^, we also 
fix the afiine freedom in u by specifying it on an initial slice S~ of Ti, which is located at an 



early time approximating the asymptotic equilibrium of the white hole at past time infinity 
/~ . The outgoing null hypersurface J~ emanating from S~ approximates past null infinity 
X~ . In Sec. [IV B|, we discuss the nature of that approximation. 



On ?i, the affine tangent to the generators rf-da = du (see Fig. |]) satisfies the geodesic 
equation n^Vtn'^ = and the hypersurface orthogonality condition n'^'V^n'^^ = 0. Following 
the approach of Refs. ||l|,@l, we project 4-dimensional tensor fields into H using the operator 

P', =6i + nj\ (2.2) 

where la = — VqW, and we use the shorthand notation ± T^ for the projection (to the tangent 
space of Ti) of the tensor field T^. The projection ± has gauge freedom corresponding to the 
choice of affine parameter u. However, the action of _L on covariant indices is independent 
of this freedom and equals the action of the pullback operator to TC. 

In addition to the intrinsic geometry of Ti, the necessary characteristic data consist of 
the extrinsic curvature quantities given by ± VJb (with gauge freedom corresponding to 
the affine choice of u). Since VJb — V^/q = 0, the independent components are determined 
by _L V(a/b), which has the decomposition 

-L V(a4) = -la^b - ^Jb + Tab- (2.3) 

Here Tab describes the shear and expansion of the outgoing rays and satisfies n"-Tab = on 
7i. Following Hayward |jl2[, we call Ua =-L nJ'Vbla the twist of the affine foliation of 7i. 



The twist also satisfies n°'uja = 0. Unlike Tab, the twist is an invariantly defined extrinsic 
curvature property of the u = const cross sections of TC, independent of the boost freedom 
in the extensions of n^ and la subject to the normalization n°'la = —1. 

Note that it is natural geometrically to associate extrinsic curvature properties of 7i with 
(_L Va)n^, in analogy with the Weingarten map for a non-degenerate hypersurface W^. The 
normalization n°'la = — 1 then leads, via Eq. (p.3|), to 



/fe(± Va)n' = -UOa. (2.4) 

Thus the twist describes an extrinsic curvature property associated with n°- as well as la- 
The other components of (_L 'Wa)'n^ are determined by the shear and expansion of Ti (which 
vanish in the special case of an isolated horizon). 

We consider the double null initial value problem based upon data on the horizon Ti 
and the outgoing null hypersurface J~ emanating from S^ , with the evolution proceed- 
ing along the outgoing null hypersurfaces Ju emanating from the u = const foliation Su- 
In this problem, the complete (and unconstrained) characteristic data on Ti are its affine 
parametrization u and the conformal part Iiab of its degenerate intrinsic metric. Similarly, 
the characteristic data on J'~ are its affine parametrization A and its intrinsic conformal 
metric Iiab- The remaining data consist of the intrinsic metric and extrinsic curvature of 



S^ (subject to consistency with the characteristic data) |jT0|,|T3[. In Sachs coordinates, this 
data on S~ consists of r, f, Ua and r_A (which determines the expansion of the outgoing 
null rays) on S~ . That completes the data necessary to evolve the exterior spacetime. In 
carrying out such an evolution computationally, the first step is to propagate the data given 
on S~ along the generators of H so that it can be supplied as boundary data to the exterior 



evolution code. This first step is accomplished by means of certain components of Einstein's 
equations. 

Einstein's equation decompose into (i) hypersurface equations intrinsic to the null hyper- 
surfaces Jui which determine auxiliary metric quantities in terms of the conformal metric 
Hab'-, (ii) evolution equations determining the rate of change dufiAB of the conformal metric 
of J'u, and (iii) propagation equations which are constraints that need only be satisfied on 
Ti. One of the propagation equations is the ingoing Raychaudhuri equation Ruu = which 
determines the surface area variable r along the generators of 7i in terms of initial conditions 
on S~ according to 

f = ^h^^'hAB. (2.5) 



The value of r on S~ measures the convergence of its ingoing null rays and Eq. ( p.5[ ) implies 
f < 0. The remaining propagation equations Rau = propagate the twist Ua along the 
generators of TC. Our coordinate conditions imply that Uu = ujx = and 

u^A = -IdxigAsW). (2.6) 

The non-vanishing components propagate according to 

{r'uAj = t'DaC-) - h^'^DBir'hAc), (2.7) 

r 2 

where Da is the covariant derivative associated with Kab- Here h'^^hBc = ^c ^^^ h^^liAB = 
0. Having determined r, this equation can easily be integrated to determine ua on 7i in terms 
of initial conditions on S~ . Once the propagation equations are solved on 7i to determine r 
and Ua, the Bianchi identities ensure that they will be satisfied in the exterior spacetime as 
a result of the J'u hypersurface equations and the evolution equations. 

The propagation of Tab (the outward shear expansion and shear) along H requires the 
Rab = components of Einstein equations, given by 



Rab = hAB y— - D^Da logrj - Daujb - Dbuja - ^ujaujb 
'-{ujAdBT + ujBdAr - hABh'^^^cdnr) 



r 

+ r'^dxfiAB - -^h'^^ih-AcdxhBD + hBudxhAc) 
+ 2{rdxr + rdxr)hAB + r{dxr)hAB + rrdxhAB- (2.8) 

where TZ represents the curvature scalar of the metric Hab- On H, the equation Rab = 
decomposes into the trace 

dudxir^) = +D^Da logr - ^7^ + D^ua + h^'^UA^JB (2.9) 

and, by introducing the dyad decomposition Hab = ''^(a^b), the trace-free part 
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m^m^ (rdxd^irhAB) - 2Daujb - 2ujaujb + -ujaDbA = 0. (2.10) 

The trace equation propagates the outgoing expansion of the fohation iS„, as determined 
by r A- The trace-free part is the evolution equation for the data Has on the fohation jZu, 
apphed at H to evolve hAB,x, which describes the shear of the outgoing rays. The outward 
shear constitutes part of the extrinsic curvature of S~ . Although its value is implicitly 
determined by the null data Has on J'^ , we view it as part of the initial data that must be 
specified on S~ . 

In summary, the data for the double null problem includes the conformal metric on Ti 
and the quantities r, r, Ua, r^x and /iab,a on S~ . Equations ( pTSf ), ( ^TTf ), ( p. 91) and ( |2.10|) are 



then used to propagate the data on S~ to all of Ti. The choice of foliation of the horizon is 
an important but complicated aspect of this problem. As in the case of Cauchy evolution, 
gauge freedom in the foliation introduces arbitrariness into the dynamical description of any 
black hole process and, in particular, the pair-of-pants structure underlying a black hole 
merger. In the case of the horizon, the natural choice of an affine foliation u removes any 
time dependence from this gauge freedom but there remains the affine freedom u —>■ Au + B, 



where A and B are ray dependent. In Sec. |3, we use the asymptotic equilibrium of the 



white hole as m — > — oo to fix this freedom. We review in Sec. |T| how the data Hab and r 
on Ti describing a white hole fission (binary black hole merger) are provided by a conformal 
horizon model [|1],0]. The additional data required on S~ can be inferred from the asymptotic 



properties of the white hole equihbrium at /", as discussed in Sec. |^. The evolution of 
the exterior spacetime by the PITT code requires a transformation of the data on 7i to a 
Bondi-Sachs coordinate system, as described next in Sees. [II B| - |1ID| . 

B. Bondi-Sachs coordinates 

The transformation to the Bondi-Sachs coordinates x°' = {u, r, x^) used in the PITT 
code consists in substituting the surface area coordinate r for the affine parameter A. Since 
the horizon Ti does not in general have constant r it does not lie precisely on radial grid 
points. Consequently, the assignment of horizon boundary values must be done on the grid 
points nearest to the boundary. Thus an accurate prescription of boundary conditions in 
the r-grid requires a Taylor expansion of the horizon data. 

In Bondi-Sachs variables, the resulting metric takes the form 

rfs^ = - ( e^^- - r^hAB^^uA du^ - 2e^^dudr - 2r'^hABU^dudx^ + PhABdx^dx^. 

(2.11) 

In relating this to the Sachs metric Eq. (|2.1| ), it is simplest to consider the contravariant form 
in which only the (?'"" components differ. In terms of the corresponding metric variables, 

2r,r„ + ^/.^^ (2.12) 

(2.13) 
(2.14) 



grr 


= e-'^^ ={r,xyW-2r^xr,AW^ 


grA 




gru 


= -e-2^ =-rA. 



Restricted to 7i, our choice of lapse and shift imply 



/3 = -^lnr;, (2.15) 



e 



2/3 



U'' = ^r>^^ (2.16) 



J.2 ' 



e 



2/3 



V = -2 rr^u + — r ^r b/i^^. (2.17) 

r 



The structure of the J7u-hypersurface equations and the evolution equations ||14| , p!5[| re- 
veals the horizon boundary data necessary for characteristic evolution. The hypersurface 
equations are 

Ar- = ^rh''''h''''hAB,rhcD,r. (2.18) 

Id 

{T^e-^PhABU^r),r = 2t' (r-^/?,^)^^ - T^h^'' DchAB,r (2.19) 

2e-^^Vr = n- 2D^DAf3 - 2D'^(3Da(3 + r-^e-^^DA{r^U^),r 

- K^e-'^hAsUp^, , (2.20) 

where here Da is the covariant derivative and TZ the curvature scalar of the 2-nietric Hab 
of the r = const surfaces (which differ from the corresponding quantities on the A = const 
surfaces). These equations form a hierarchy which can be integrated radially in order to 
determine /3, [/^ and l^ on a hypersurface Ju in terms of integration constants on Su, once 
the null data Hab has been evolved to Ju- 

The evolution variable Hab can be recast as a single complex function, since det(/?,AB) = 
det{qAB) = q{x^) is independent of u and r. The code treats such functions on the sphere 
in terms of stereographic angular coordinates based upon the auxiliary unit sphere metric 
qAB- Tensor fields are represented by spin- weighted functions using a computational version 
of the 9-formalism [|I6| based upon a complex dyad g^, satisfying qAB = <l(AqB)- (Note 
that this departs from other conventions |T^ in order to avoid unnecessary factors of \/2 
which would be awkward in numerical work.) For example, the vector field va is represented 
by the spin- weight 1 function v = q^VA- Derivatives of tensor fields are represented by 5 
operators on spin-weighted functions by taking dyad components of covariant derivatives 
with respect to the unit sphere metric. Our conventions are fixed in the case of a scalar field 
\E' by S\E' = q^V A^ ■ In these conventions, 

(8a - m)r] = 2sr] (2.21) 

for a spin- weight s field rj. 

The conformal metric Hab is represented by the dyad components J = hABq^q^ /2 and 
K = hABq^q^ /2, with the determinant condition implying K"^ = 1 + J J. As discussed 
in W^, the curvature scalar corresponding to Hab is 



n = 2K- mx + i[8V + g2 J] + ^[gjgj - SjgJ]. (2.22) 

2 4i\ 

In terms of J, the evolution equation takes the form 
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2(rJU-(r-V(rJ),)^^ = ^, (2.23) 

where Sj is explicitly given in term of spin- weighted fields in Ref. [Q. Sj contains only first 
derivatives of J and predetermined hypersurface quantities so that it does not play a major 
role in the integration of Eq. (|2.23| ) on a given null hypersurface J'u- 

Similarly, in integrating Eq. ( p.l9|) for U"^, the code uses the spin-weight 1 fields U = 



qAU^ and Q = q'^Qa, where 

Qa = r^e-^^hABU^, (2.24) 

so that Eq. ( p.l9|) reduces to the first order radial equations 

{r'Q)^^ = q^[2r' (r-^A^) , - r'h''^ DchAsA. (2.25) 

U,r = r-'e'^qAh^''QB, (2.26) 

with the right-hand sides then rewritten in the S-formalism in terms of spin-weighted fields. 

The integration constants necessary to evolve Eqs. ( p.l8| ) - (|2.23| ) are the boundary 



values of r, J, /?, Q, U and V on Ti.. Part of these boundary conditions are determined 
by the characteristic data on Ti and the remainder are determined by gauge conditions and 
the solution of the propagation equations. When the horizon Ti has constant r (as in the 
Schwarzschild case), this is precisely the data necessary to initiate the radial integrations 
in the PITT code. In the generic case, Ti does not lie on grid points and the initialization 
of the r-grid boundary also requires drJ, dr(3 and drV on Ti to provide a Taylor expansion 
consistent with a second order accurate code.. 



C. From horizon variables to Bondi metric functions 

We now describe the calculation of the Bondi-Sachs metric variables in a neighborhood 
of the horizon as necessary for a characteristic evolution of the exterior space-time. The 
essential piece of free horizon data is the conformal metric Hab or equivalently the metric 



function J. In Sec. |T|, this data is supplied for the case of a binary horizon by the conformal 
model, in which case the value of r on the horizon is also supplied. In a more general 
situation, where only J is prescribed, r can be determined from its value and time derivative 
r on S~ using the ingoing Raychaudhuri equation ( |2.5D , which in spin- weighted form is 

r = -'-{JJ-K'). (2.27) 

Here, we have used the dyad expansion of the conformal metric h"^^ which can be written 
in terms of J and K as 

2 h^^ = -Jq^q^ - Jq^f + K [q^f + q^q^) . (2.28) 

Once both the intrinsic geometry J and the radial coordinate r are known, the metric 
functions /3, U and W (as well as their radial derivatives) can be calculated using a hierarchy 



of horizon propagation equations [similar to the hierarchy of hypersurface and evolution 
equations (p.l5|) - (|2.23|) used to propagate fields along a Ju null hypersurface]. 



Next in this hierarchy is the propagation Eq. ( p.7| ) for the twist uj = q^oJA , which has 
spin-weighted form 

(r^iu)- = r^Q{-) + -[Jd{r^k) + Jd{r^j) - K^{r^j) - Kd{r''k)] 

2 

+ ^{JQJ - 4:kQK + 2kdJ + 3 JSJ - 2JQK) (2.29) 

and determines u given its initial value on the slice 5*" and the conformal horizon geometry 
encoded in J. 

Next, the value of f3 follows from Eq. ( |2.15D with r ^ determined from Eq. ( |2.9| ) in the 
spin-weighted form 

dudxir^) = -g[i^(Slogr + u)] + -Q[K{B logr + u)] - -8[J(8 logr + u)] - -d[J{Q log r + u)] 
- -n--Juj^--Juj^ + Kujiu, (2.30) 

where Eq. (|2.22| ) supplies the spin-weighted form of the curvature scalar TZ of the metric 
Hab- This determines r a on the entire horizon given its value on the initial slice S~ . 

Next, the Bondi metric functions U and V are evaluated on the horizon using Eqs. (p. 16 ) 
and ( |2.17| ) in the spin-weighted forms 



e2/3 
f/ = ^- [jQr - KQr) (2.31) 

2/3 

V = -2rr^u + ^ (-J i^rf - J {QrY + 2 K (Qr) (Br)) . (2.32) 

In summary, the Bondi metric functions r, J, (3, U and V are determined by the hierarchy 
of Eqs. (|2.27|) and ( p.29| ) - (|2.32|) . When the value of r is given, as in the conformal model, 



Eq. ( p.27|) is not needed. 

D. Extending the Bondi metric off the horizon 

The evaluation of the Bondi metric functions on radial grid points in the neighborhood 
of the horizon requires the r-derivatives of J, /3, U and V. For this purpose we note that, 
on any given ray, the quantity drF is known once both dxF and r^\ are determined, e.g. 

Or J = dxJ/r^x = e^^dxJ. (2.33) 

We assume below that r^x has already been determined on the horizon by integrating 
Eq. ( p^Ol) . 



We obtain J,, on the horizon (and hence K^^) in terms of dxJ, which is determined from 
its initial value on S~ by integrating Eq. ( |2.10| ) in the spin-weighted form 
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= 2r'^dudxJ + 2rrdxJ + 2rjdxr - Jr^JdxJ + JdxJ - 2kdxK) 

- {1 + K^){duj + Lo^ ~ 2ujQ log r) + uj{JdK - KdJ) + uj{KQJ - JQK) 

- j{jm - K{m + Scu) + Juj"^ - 2Kcucu + 2K{ujElogr + coQhgr) - 2JtuBlogr). (2.34) 



We calculate P^r from the Raychaudhuri equation for the outgoing null geodesies, Rxx 
0, which in Sachs coordinates takes the form 



^ACuBDi 



r,xx = --h'''-h''''hAB,xhcD,x 



with spin-weighted version 



r,AA = ^ {{K^xf - J,xJ,x) 



(2.35) 



(2.36) 



which, together with Eq. (|2.15| ), gives 



M T o V ^ ^ ^^ T ) ' 



(2.37) 



We obtain U^r (or Q) from the twist, which in Bondi coordinates takes the form 

UA = dAl3 + -dAT + ]-dr{hAB)h''''dcr - "^e-^^ hAsdrU'' , 
r 2 2 

with spin-weighted version 



(2.38) 



uj = Qf3 + -Qr + hdrJ)iK(5r - Jdr) + ^{drK){KQr - Jdr) - ^. (2.39) 

This determines the boundary data for Q in terms of quantities that are already known on 
the horizon. Given Q, the value of drll follows from Eq. ( |2.26| ), which has spin-weighted 
form 



,2/3 



t/., 



— {KQ - JQ) . 



(2.40) 



Finally, we compute drV by obtaining Vx from the A-derivative of Eq. ( p. 121) , 

p-2/3 



V,-V{2p,x + ^-f 



ir,x)lW + {r,xfW,x - 2 (r Ar,^) ^ 



W^-2r.r,Wi 



O o , o^^Ar,B,AB c.r,xr^Ar\B,AB 
2rxxr,u-2rxrxu + 2 — h -2 h 



, r\Ar\B.AB 



(2.41) 



The first three terms in the right-hand side vanish on 7i due to the gauge conditions. We 
express the others in spin-weighted form using Eq. ( |2.6| ) to obtain 



-^^^AWJi=-2 



(— JcjSr — Judr + K (cjSr + tuSr)) 
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(2.42) 



and Eq. ( p.2^ ) to obtain 

2r A^r b/i^^ = -J{dr^x){^r) - J(8r,A)(Br) + K ((5r,A)(8r) + (Br,A)(5r)) (2.43) 

2r ^r b/i1^ = -Jxi^rf - Jx{^rf + 2 K^x{Qr){dr). (2.44) 



Then Eq. ([21411 ) gives 



/ e"2^\ 2 - - 

V^ = V i2p^x + ) +4rr„/3,A-2rr,A„ + - {-Jcodr - JcuQr + K {cudr + coQr)) 

+ — (-J (Sr a) (Sr) - J (Sr a) (Sr) + iT (gr,A) (Sr) + K (Sr a) (9r)) 

+ \ (j {Qrf + J {drY - 2 K (dr) (8r)) 

+ — (-Ja (Sr)2 - Jx (Br)' + 2ir,A (5r) (8r) j (2.45) 

in terms of previously determined quantities on the right-hand side. 

The values of each metric function (J, f3, U, V) and its first radial derivative can then be 
used to consistently and accurately initialize r-grid points near the horizon. This in turn 
allows the evolution code to determine the entire region extending from the horizon to X^, 
as long as the coordinate system remains well behaved. 

III. CONFORMAL MODEL FOR THE AXISYMMETRIC HEAD-ON COLLISION 



The conformal horizon model |l|J^ supplies the conformal metric h^B constituting the 
null data for a binary black hole. The conformal model is based upon the flat space null 
hypersurface 7i emanating normal from a convex, topological sphere Sq embedded at a 
constant inertial time t = in Minkowski space. Traced back into the past, 7i expands 
to an asymptotically spherical shape. Traced into the future, Ti pinches off where its null 
rays cross, at points X, or where neighboring null rays focus, at caustic points C. Figure |^ 
schematically illustrates the embedding of 7i in Minkowski space for the case when Sq is a 
prolate spheroid. The spheoridal case generates the horizon for the axisymmetric head-on 
collision of two black holes, or here in the corresponding time reversed scenario, the horizon 
for an initially Schwarzschild white hole which undergoes an axisymmetric fission into two 
white hole components. The white hole horizon shares the same submanifold 7i and the 
same (degenerate) confomal metric as its Minkowski space counterpart but its surface area 
and afiine parametrization differ. As a white hole horizon, 7i extends infinitely far to the 
past of Sq to an asymptotic equilibrium with finite surface area. 

The intrinsic white hole geometry of Ti is obtained by the following construction based 
upon the fiat space null hypersurface. A t foliation of 7i is induced by the t foliation of 
Minkowski space, with t and the Euclidean coordinates (x, y, z) determining a Minkowski 
coordinate system. The prolate spheroid Sq is given by 

x'' + y^ + - = a^ (3.1) 
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with e > 0, or alternatively in angular coordinates y"^ = {rj, 0) by 

X = asinr] cos (3.2) 

y = a sin r] sin (3.3) 

z = y/l + eacosT] . (3.4) 

In these coordinates, the Minkowski metric induces on Sq the intrinsic metric 

QABdy^dy^ = a^ ( (l + e sin^ r^) drf + sin^ i] d(j)^) . (3.5) 

The principal curvature directions on So are the polar and azimuthal directions, with 
corresponding radii of curvature 

^ = a 7r^= (3-6) 

and 

(l + esin^r^)!/^ 

r<i^ = a j=== . (3.7) 

Vl + e 

Each generator of H encounters two caustics at t = r^ and t = r^ (or one degenerate caustic) 
if continued into the future but, along a typical generator, 7i first pinches off at a cross over 
with another generator before a caustic is reached. 

The time dependent metric of the t foliation is given by 

gABdy^dy^ = a^{(l- — \ {I + esin^ ri)drf + (l - — \ srn^ r]d<pA . (3.8) 

As t ^ — oo, g^^/grfq -^ sin^r7[l + esin^?7] so that the conformal metric in these coordinates 
does not approach the standard form of the unit sphere conformal metric. For this purpose, 
it is convenient to introduce new angular coordinates x^ = {9, 0) in which the conformal 
metric asymptotes to the unit sphere metric qAsdx'^dx^ = dO"^ + sin^ 6*^0^. This requires 
that 

^^ - ^^ (3.9) 



sine' sinr^Y 1 + esin^r^ 
with the solution 



tan 6* = Vl + etanr]. (3.10) 

(Here the boost freedom corresponding to the unit sphere conformal group has been fixed 
by requiring the transformation to have reflection symmetry about the equator.) 

In these a;^ coordinates, gAsdx^dx^ = f'^hAndx^dx^ , where det{hAB) = det(gAB) = q, 

(l + ecos^a*)^ V ^0/ V ^<^/ 
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and 



^^-"(l + ecos2 0)3/2 (^-^2) 



:i + ecos2 0)1/2- ^3.13) 



In the prolate case (e > 0), r^ > r^. 

We apply the conformal horizon model to endow Ti with the intrinsic metric qab = ^"^(Jab 
of a white hole horizon. The conformal factor Q is designed to stop the expansion of the 
white hole in the past so that the surface area asymptotically hovers at a fixed radius. As a 
result, hAB{t-i x"^) = hAB{t-, x^) is the intrinsic conformal metric of the horizon (as well as of 
the flat space null hypersurface) . With the dyad choice q^ = (1,2/ sin 0), Eq. ( p.8|) implies 
that the spin-weight-2 field 

J = i/,«ft.B(<. .-) = i (4^) - i (l^) (3.14) 

2 2 \t — r^^ 2 \t — rg^ 

is the conformal null data for the white hole horizon in the t foliation. 

The surface area of the white hole is related to the corresponding surface area of the fiat 
space null hypersurface hj r = Qf. A conformal factor with all required behavior to produce 
a non-singular white hole is given by |l|J^ 

2 

n = -R^iu + —^—)-\ (3.15) 

I2[p — u) 

where R^o is the initial equilibrium radius, p is a model parameter (designated by p in 
Refs. |l|,@]), 0" is the difference between the principal curvature radii, 

a\e\ sin^ 9 

and u is an afiine parameter along the generators of Ti with the same scale as t but with 
origin u = chosen to lie midway between the caustics, i.e. u = t — tq where 

{re + r^) 2 + 2e-esin20 

2 2(1 + ecos^ ayi'^ 

is the mean curvature of Sq. For an initially Schwarzschild white hole, i?oo = 2M. 

As a 3-manifold with boundary, Ti represents both the white hole horizon and the flat 
space null hypersurface. Both extend infinitely into the past and continue into the future 
to the boundary where Ti pinches off at caustics and cross over points. Smoothness of the 
white hole requires that the parameter p > ctm/vTs, where cta/ is the maximum value of a 
attained on Sq. For a prolate spheroid, the maximum occurs at the equator and o"m = a\e\. 

TC must obey the Raychaudhuri equation ( p3| ) both as a fiat space null hypersurface and 
as a white hole horizon. Both geometries have the same intrinsic conformal metric so that 
they must have the same rate of focusing in terms of their respective affine parameters, in 
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accord with the Sachs optical equations ||T8[ . As a result of the different behavior of their 
surface areas due the conformal factor Q, the affine parameter t on the white hole horizon 
is related to its flat space counterpart t according to |1| 

di (12M(n-p)-a2)2(5j9-;U-2u)2(2p/M-i) ' ^'^- ^> 

where 

fx = Vl3p2 - a2. (3.19) 



Here the affine scale of t is fixed by the condition dt/dt ^ 1 as t ^ — oo. Equation ( p.l8| ) 
determines the deviation of a slicing adapted to an affine parameter of the white hole horizon 
from the original slicing given by the Minkowski embedding. The angular dependence of 
the crossover time t^ at which the horizon pinches off leads to the change in topology of the 
white hole associated with a pair-of-pants shaped horizon. From Eq. (|3.13|) , the pinch-off 
occurs at tx = r^ = a/ a/1 + ecos^ 9, or 

aesin^ 9 . . 

ux = -777^ 2m^- 3.20 

2(1 + ecos"' 9y^''^ 

Consequently, ux as well as a vanish at both poles so that integration of Eq. (|3.18| ) implies 
that t;^:' — > cxD at both poles. This creates the two legs of the pair-of-pants since t^ remains 
finite at all other angles. 

Note that A = A{u/a,p/a, e) within our conformal model, so that t = t(t/a,p/a, e). The 
same scale dependence on a also holds for the horizon variables f and J, as is evident from 
Eqs. ( p.ll|) and (|3.14|) . In addition Eq. ( p.l5|) implies r = r{Roo/a, u/a,p/a, e). Thus we can 
scale a = 1 without any loss of generality of the model. 

IV. THE CLOSE APPROXIMATION AND DATA ON 5 

In addition to the conformal metric on 7i (discussed in the last section), specification of 
the radius and extrinsic curvature of the initial slice S~ completes the necessary data on Ti. 
Our strategy is to locate S~ at an early quasi-stationary era and approximate these data 
by their equilibrium values as a Schwarzschild white hole. In the linearized approximation, 
the conformal metric of 7i corresponds to a perturbation of the Schwarzschild background. 
However, a comparison of how the fully nonlinear data for the head-on collision deviate from 
their linearized counterparts reveals features which can not be described perturbatively. The 
chief geometrical issues to be discussed here are the asymptotic properties of the horizon at 
I~ , the behavior where it pinches off and the analogue of the bifurcation sphere occurring 
in a Schwarzschild horizon. As indicated in Fig. |^, we choose S" to correspond to an early 
Minkowski time t = t_ so that it is initially quasi-spherical for all eccentricities e of the 
ellipsoid Sq located at t = 0. The criterion that the initialization be at an early time is 
|t_| >> vq. In that approximation dt/dt ^ 1 on S~ and Eq. ( |3.14| ) gives 



J^ ^ '-^t^^. (4.1) 

t_ 
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A. Perturbing a Sch^varzschild horizon and the close approximation 

In the Cauchy treatment of the close approximation to the axisymmetric head-on colhsion 
of black holes , the background spacetime is the exterior Kruskal quadrant of the extended 
Schwarzschild space-time. The trousers shape of the binary horizon is beyond the scope of 
a perturbative treatment. Here we present a fully nonlinear characteristic treatment of the 
(time reversed) head-on collision as a 1-parameter sequence of binary collisions with the 
Schwarzschild case as a limit. A trousers-shaped horizon exists for each non-Schwarzschild 
member of the sequence so that it is possible to investigate its behavior in the characteristic 
version of the close approximation. 

The conformal horizon model of a Schwarzschild horizon is the highly degenerate case in 
which So is geometrically a sphere. Then hat = Qab, J = and r = _Roo = 2M, so that the 
horizon is stationary. In the {u, A, x^) coordinates of the Sachs metric, the Schwarzschild 
geometry is described by Eq. ( |2.1[ ) with W^ = 0, r = 2M — Am/(4M) and 

W = . (4.2) 

AM-8M2 ^ ^ 

These coordinates cover the entire Kruskal manifold with remarkably simple analytic be- 



havior, as observed by Israel [19|. On the white hole horizon, given by A = 0, m is an affine 
parameter with its origin fixed so that m = on the r = 2M bifurcation sphere where 
dxr = 0. 



The spin-weighted versions ( |2.27| ) - (|2.34|) of the Einstein Eqs. (|2.5| ) - (|2.1(]| ) lead to the 



following linearized equations governing the perturbation of a Schwarzschild horizon. The 
ingoing Raychaudhuri Eq. ( p. 27 ) simplifies to 



f = 0, (4.3) 

so that we can set r = 2M on 7i, where M is the background Schwarzschild mass. Equa- 
tion ( p.29|) then reduces to 



so that 



(r^oo}=-lr^dj, (4.4) 



u = -SJ/2, (4.5) 



where we fix the constant of integration by requiring that the perturbed black hole come to 
equilibrium as a Schwarzschild black hole with J = to' = OasM— i> — cxd. Next Eqs. ( p.30| ) 
and (|2.34| ) reduce to 



2 4 
and 



rd^dxr = -1 - 1(82 J + g" J) (4.6) 



= 2rdudx{rJ) + SSJ + J. (4.7) 
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It is convenient to set r = prM-, where 

TM = 2M - ^ (4.8) 

is the background Schwarzschild value and where, restricted to H, p = 1 + O(J^). Then, to 
first order in the perturbation, Eq. ([4 .61) reduces to 



and Eq. ( [4.7[ ) reduces to 



r'9„9Ap = -^(eV + g2J) (4.9) 



Q = 2r'^dudxJ + ^^J -uj. (4.10) 

Note that in the perturbative hmit the right-hand sides of of the equations for the extrinsic 
quantities u;, r ^ and J\ only depend on the intrinsic geometry. 

As the parameter e describing the eccentricity of the spheroid approaches zero, the binary 
black hole horizon approaches the Schwarzschild horizon of a single spherically symmetric 
black hole. Thus the close approximation for a binary black hole can be described by a 
perturbation expansion in e. In order to integrate the perturbation Eqs. ( [4.9| ) and ( [4. 101) , 



with J supplied by the conformal model, we must relate the affine parameter t of the 
conformal model to the restriction to Ti of the coordinate u of the Sachs metric. We set 
t = Au, where 

A= n^ \n,3/. ='^ + 0ie) (4.11) 

(1 + ecos^ tly/'^ 

fixes the relative scale freedom in the affine parameters so that the perturbation has the 
early time behavior 

aesin^^ ,, ^^. 

J (4.12) 

u 

in accord with the quadrupole nature of the close approximation. 



B. Data on S 

The conformal model described in Sec. |T| supplies the values of J and r on the entire 
horizon. Other quantities have to be initiated at an early cross section S~ , near which the 
horizon behaves as a perturbed Schwarzschild horizon. As indicated in Fig. |^, we locate 
S~ at a constant Minkowski time t = t_, as well as at a constant white hole affine time 
u = U-. The relationship t{u) is then determined by integrating du/dt = 1/{AA) with initial 
conditions determined on S~ and A and A given in Eqs. ( ^.1^ ) and ( [4.111) . This determines 
the horizon data J{u, x^) from J(t, x^) given in Eq. ( p.l4| ). Similarly, r{u, x^) is determined 
from r(t,x^) = fl(t,x^)r{t,x^), with f and fl given by Eqs. ( p.ll| ) and ( p.l5| ), respectively. 
The requirement that S^ be located in the quasi-equilibrium era implies that r_ ^ 2M to 
an excellent approximation. 
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The way in which the values of U- and t_ determine the location of So relative to S~ 
depends upon the parameters entering the conformal model. Even in the limit e ^ 0, where 
the conformal model yields a Schwarzschild horizon, this relation depends upon the model 
parameter p. The simplest limiting case is when p and e both vanish. Then A = A = 1 so 
that t — t^ = u — u^ and the Minkowski spheroid So (where t = 0) is located dX u = U-—t-. 
Thus a negative value of t_ locates So to the future of S~ . 

In the Schwarzschild limit, dxr -^ OxTm = —u/AM, where m = on the r = 2M 
bifurcation sphere B. This allows us to use the initial value of d\r to determine the location 
of S~ on the horizon by setting 

avr- = -^ (4.13) 

on S~ . Thus specification of the initial outward expansion oiS^ fixes the translation freedom 
in the affine parameter u. The requirement |u_| >> M ensures that the initialization be at 
an early time. 

A 1-parameter family of horizon data for a head-on collision results from choosing So to 
be an e-family of spheroids (with e > 0). The close approximation, corresponding to the 
behavior linear in e, provides insight into the asymptotic structure of the horizon at early 
times. In this linear approximation, Eq. ( |3.14D implies 



ea sin^ 9 / . , . s 

J ^ (4.14) 

t — a 

on the horizon; and Eqs. ( ^.l(j| ) and ( |3.17| ) imply 

a^easiii^e (4.15) 

and 

r^^!l±Il = a(^l + --ecos^ey (4.16) 

In addition, r = 2M + O(e^). (We use ~ to denote approximations valid for small e and ~ 
to denote asymptotic approximations at early times.) 

The early time asymptotic behavior of all the horizon variables can be explicitly evaluated 
in this approximation. Equation (|4.5|) determines the asymptotic dependence 

cj~-8J/2, (4.17) 

on the assumption of an initially non-spinning Schwarzschild white hole. Using the identities 
g2cos2 6' = 2sin^6', {d^d'^ + QW) cos^ 9 = AS^cos'^ 9 - i), the commutation relation (g^ 
and the property BSV^m = — ^(^ + l)Yem, Eq. (^^) reduces to 

At early times, where A ~ 1, this integrates to give 



r^dxp r~^Qea(cos^e-^] log (L-^\ (4.19) 



where we set the integration constant so that dxp = on 5 in accord with Eq. ( [4.13[ ). 
Similarly, at early times, the integral of Eq. ( |4.1CI| ) gives 



^ ^ 3easin^0, ft — a\ / , „„x 

*^~-^;^'os(|-3^). (4.20) 

Note that d\p and d\J have logarithmic asymptotic behavior at /~. In the exterior 
evolution code this singular behavior is renormalized by dealing instead with the quantities 
drP = d\p/d\r and drJ = d\J/d\r. Since d\r ~ —u/AM at J~, these quantities both go to 
zero as log u/m as u —>■ — oo. This justifies initializing these quantities to zero on S~ . The 
initialization error is 0{logu_/u^) and converges to as u_ ^ — oo. In summary, given 
J and r on the horizon via the conformal model, the remaining data necessary on S~ are 
initialized according to 

oo. = -^BJ., (4.21) 

dxp. = 0, (4.22) 

dxJ- = 0. (4.23) 

With this initialization, all Bondi-Sachs start-up variables for the exterior evolution are 
asymptotically well defined at /~ except for e"^'^ = d\r ~ — 'u/(4M). This is handled by the 
renormalization f3 = (Sm + I3r where e~'^^'^'' = —u/{4M) and /?/?—> at /~. In the evolution 
code, the singular part Pm is analytically factored out of the Bondi equations. The regular 
part is given by 

e-.» . 1 _ !^1!*P. (4.24) 

u 



Then, referring to Eq. ( [4.1^ ), f3p> can be initialized to on 5 with O (log «_/«_) accuracy. 



The early time approximation breaks down before reaching the crossover points X where 
t ~ a, as evident from Eqs. ( 4.18|) and ( [4.2C ). The way that the horizon pinches off at X in 



a sequence of models as e — »■ is sensitive to the behavior of the model parameter p along 
the sequence. We consider here the case p = const (independent of e as well as angle). 

The crossover points occur at Minkowski time ix = a/Vl + ecos^ 9 ~ a[l — (e/2) cos^ 9]. 
The corresponding values ux = t^/A are found from integrating Eq. ( |3.18[ ). The rays on the 
poles of the prolate spheroid, where a = 0, do not focus. Since the horizon surface area is 
asymptotically constant, this lack of focusing implies that the horizon persists forever along 
the poles, i.e. for — oo < u < +oo (independent of the value of e). 

In order to investigate the location of X along the non-polar rays, consider the small e 
behavior 

1_ (3p^ - bpii + -»^)V(5 + VT3)p - 2m \ 4/^13 

a"^ u^{u-py \{5-VTS)p~2u) ' ^' ^ 
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where 



^ -.2n\ _ X X «e . 2, 



u^t-a[l + --ecos^e] ^ t - tx - — sin^^. (4.26) 



The dominant behavior near X is revealed by the further approximation 

l_,9pV5 + v^y"^_ (4,27) 



A m2 V5_ v^/ 
which is vahd for u << p, e.g. near A" where u ^ 0. In this approximation, 

ux-u_= TT^V ^ ^^r7{ + -- ' 4.28 

so that ma" — ^ C)0 as e — > along all rays. In this limit, the entire cross over seam on the 
pair-of-pants is mapped to t = oo. As a result, the corresponding first order perturbation 
theory for this version of the close approximation is well behaved on the entire white hole 
horizon, not just the segment bordering the exterior space-time. 

Of special physical importance is the location of the crossover surface X where the 
horizon pinches off relative to the surface B defined by dxr = 0. B represents a boundary 
for the Bondi evolution resulting from the breakdown of the surface area coordinate r. In 
the e = Schwarzschild limit, B is the r = 2M bifurcation sphere located at ti = and X 
lies at M = oo (at /"'"). In this limit, the white hole fission takes place in the infinite future. 
For small but nonzero e, the dominant O(l/esin^0) dependence in Eq. (|4.28| ), implies that 
X lies at a finite but large value, with the point at the equator where the white hole fissions 
(the crotch in the pair-of-pants picture) located at the earliest point on X. However, B 
remains within 0(e) of its Schwarzschild location at m = 0. Thus, for small e the fission is 
"hidden" beyond B in the sense that it is not visible to observers at X"*". From the time 
reversed view of a black hole merger, the individual black holes would merge inside a white 
hole horizon. 



V. NONLINEAR HORIZON DATA 

The close approximation results just described for a white hole fission, when reinterpreted 
in the time reversed sense of a black hole merger, imply that the individual black holes 
merge inside a white hole horizon corresponding to the marginally anti-trapped branch of 
the r = 2M Schwarzschild surface. Of prime importance in the non-perturbative regime is 
whether an entirely different scenario is possible in which the individual black holes form 
and merge without the existence of a marginally anti-trapped surface (MATS) on the event 
horizon. The ingoing null hypersurface which intersects the horizon in such a MATS has an 
extremum in its surface area. As a result, the Bondi surface area coordinate based upon an 
ingoing null foliation is singular at the MATS. A Bondi evolution carried out backward in 
time on these ingoing null hypersurfaces would terminate at the MATS. In particular, the 
absence of a MATS to the future of the merger is a necessary condition for a Bondi evolution 
backward in time throughout the entire post-merger period. (See the discussion below for 
more technical details.) 
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It is possible to avoid this problem by means of a null evolution using an afiine parameter 
rather than a surface area parameter as radial coordinate. Null evolution codes in different 
gauges than the Bondi gauge have been developed |^^. However, at present such codes 
have not been successful in the stable evolution of dynamic horizons. 

Restated from the alternative viewpoint of a white hole fission, as being pursued here, the 
absence of a marginally trapped surface (MTS) prior to the fission is a necessary condition 
for a Bondi evolution forward in time throughout the pre-fission period. Thus, the absence of 
a pre-fission MTS is necessary in order to carry out our strategy for obtaining the complete 
post-merger wave form of a binary black hole by means of a Bondi evolution. 

The bifurcation sphere B in the Schwarzschild space-time is a MTS on the white hole 
horizon (and in this degenerate case also a MATS). As already discussed in Sec. [IV B| and 
explicitly demonstrated in Sec. |VI|, in the small e regime of our model the corresponding 
white hole fissions at a very late time well beyond the MTS. This is the expected effect of 
a non-singular perturbation: in the e — >■ limit there is only one white hole so that the 
fission is hidden at I^. The behavior in the non-linear regime is not so easy to predict. The 
following discussion of the properties of an MTS and its relation to a Bondi boundary B 



provides a basis for understanding the computational results of Sec. ^, 

It should first be noted that, unlike the definition of a MTS, the definition of B is foliation 
dependent. A MTS is a topological sphere with one non-diverging normal null direction and 
the other divergence-free. The Bondi boundary B on a. white hole horizon 7i parameterized 
by u is the earliest slice u = ub whose outward null normal is not strictly diverging at all 
points. In the relation between (m. A, x^) Sachs coordinates and (m, r, x^) Bondi coordinates, 
this implies that dxriuB^^^x"^) = at some point of B. In other fohations of ?i, B (if it 
exists) could occur later or earlier than in the afiine foliation considered here. 

In the afiine foliation, a MTS on of Ti can be described in Sachs coordinates in the form 
u + F^x"^) = 0, A = 0. The outgoing null normal to the MTS is 

La = -adaX-da{u + F) (5.1) 

where 

1 



« = ^X''^^AF)dBF (5.2) 

which, with n°'da = du, defines the projection tensor 



The MTS satisfies 



which has the coordinate form 



1 



^a6_^ab^2L(V). (5.3) 



r'VaU = (5.4) 



r'VaL, = --[D''DAF 



AB\ 



{D''F)dx9uA + aduir') - dx{r') - (9„/i^^)(D^F)Z}Bi^]U=-F 

(5.5) 
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(where DAhsc = 0). The MTS, if it exists, can be located by solving Eq. ( p.5|) . 

The following two propositions relate the Bondi coordinates to the existence of a MTS 
on a white hole horizon 7i: 

Proposition I. A Bondi cross section u = const satisfying d\r = is a MTS. 

Proposition II. A MTS cannot exist in a region u < ug in which dxr > 0. 

The first proposition follows immediately from setting F = in Eq. (|5.5|) . The second 
proposition follows from noting that at some point on the MTS the function F would have 



a maximum where D^F = and where Eq. (|5.5| ) would reduce to 



D^DaF = dx{r'^) > 0. (5.6) 

But the inequality in Eq. ( |5.6| ) precludes the existence of a maximum. 

The second proposition establishes that a MTS cannot form before the Bondi boundary. 
Thus a Bondi evolution might terminate prematurely due to an injudicious choice of foliation 
of TC. A computational module for locating a MTS on a null hypersurface has been developed 
and successfully used for long term tracking of a moving black hole |^. In future work, this 
module will be applied to binary horizons. Here we consider only the less geometrical Bondi 
boundary B. 

As e increases into the nonlinear region, the effect on B can be seen from integrating 
Eq. ( p.9|) over the sphere. In doing so, we note that the terms which are divergences integrate 
to zero and we can apply the Gauss-Bonnet theorem to the curvature scalar term to obtain 

du (f dx{r^)dS = -47r + / h^^cuAtOBdS. (5.7) 

The term —An is responsible for the formation of B in the background Schwarzschild case. 
The nonlinear correction due to the twist is of an opposite sign and delays the formation. 
Thus, as e increases from 0, the formation of B is delayed while the location of X moves to 
earlier times. Although Eq. (|5.7| ) only describes averaged angular behavior, it suggests that 
sufficient nonlinearity might cause the white hole fission, located on the equator of A", to 
occur prior to B . (In the time reversed case, this would alow the merger of individual black 
holes without the necessity of a MATS and the consequent singularity in its past implied by 
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VI. NUMERICAL RESULTS 

The scenario hypothesized at the end of Sec. |V|can indeed be demonstrated by integrating 
the equations underlying the conformal horizon model. At an early time, the equilibrium 
conditions on the white hole horizon imply that r = 2M and dxr = —u/AM > 0. As the 
horizon evolves, the surface area r decreases along all rays but, for the axisymmetric and 
refiection symmetric fission considered here, it decreases fastest along the equatorial rays 
where the pinch-off first occurs. The outward expansion measured by Qqut = 2dxr/r also 
initially decreases along all rays, although this process can be reversed by the growth of 
nonlinear terms, as indicated by the ray-averaged behavior governed by Eq. ( p.7| ). In the 
close approximation, the expansion goes to zero along all rays before the horizon pinches 
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off, i.e. the crotch at the center of the pair of pants is hidden behind a MTS. The crucial 
question in the nonhnear regime is whether the horizon can pinch off before the formation 
of a MTS, i.e. whether the crotch is bare. For the related question in terms of a Bondi 
boundary rather than a MTS, the issue is who wins the race toward 0, the radius r or the 
expansion Qqut along some ray. 

We conduct this race for each of a sequence of models in the range < e < 10~^, 
with the remaining parameters fixed at M = 100, U- = —100, t_ = —10, a = 1 and 
p = (10~^ + 10~^)/v^l3 (just above the minimum value of p allowed by regularity of the 
conformal model for this range of eccentricities). 

We monitor the minimum value over the sphere of the expansion of the outgoing null 
rays on the horizon, and of the Bondi radius of the horizon. The results are displayed in 
Fig. ^ in terms of values of r and Qqut normalized to 1 at the initial time, so that the race 
starts out even. 

Panel (a) of the figure, for the small value e = 10^^, shows little deviation from a 
Schwarzschild horizon. A Bondi boundary B forms at m ?« as a consequence of the zeroth 
order in e term 7^ ^ 2 in Eq. ( p^ ) (7?. = 2 for a unit sphere), which causes Qqut to 
decrease linearly with u. We have verified that the initial slope of Qqut as seen in the graph 
corresponds to the expected Schwarzschild value. 

For e = 10~^, still near the close limit, the radius of the horizon hardly changes before 
the Bondi horizon forms, as illustrated in panel (6). However, the deviation of the expansion 
from a pure linear-in-time behavior is noticeable and its deviation from spherical symmetry 
as a function of ray is also noticeable in the full numerical data. (The angular behavior of 
the relevant geometrical quantities is discussed more fully below.) As e increases to 10~^ 
in panel (c), both the radius and the expansion show markedly nonlinear behavior but the 
expansion still readily wins the race toward zero. However, its margin of victory gets smaller 
with increasing e as manifest in panel (d) for e = 10"'^ in which the race is nearly a tie. 

For e = 10"'^, as shown in panel (e), the radius now wins the race. The fission takes 
place while the expansion is still significantly large, at about 90% of its initial value. For a 
larger value of e, as shown in panel (/), the effect is even more dramatic. The radius makes 
a sudden plunge to zero to win the race before the expansion has undergone any appreciable 
change. This is a dramatic nonlinear effect. Although the radius and expansion begin the 
race at the same starting point (in rescaled units), the expansion begins with a flying start 
(its initial slope in the figures) and gets accelerated by linear effects whereas the radius 
starts from rest and only gets accelerated by quadratic or higher nonlinearities. 

Figures ^ and ^ are embedding pictures of the horizon which reveal the angular behavior 
of the race. (The construction of the embedding pictures is explained in Ref. |0). The 
darkened portions of the pictures indicate where the expansion has gone negative. Figure |^ 
shows that when the expansion reaches zero first it does so at the pole; whereas the radius 
first reaches zero at the equator where the horizon pinches off to form separate white holes. 
Figure ^ shows that for e = 10~^ the pinch-off occurs before the outward expansion has gone 
to zero along any ray. 

More insight into the angular behavior is provided by surface plots as functions of (m, 9) 
of the quantities r, Qqut, the curvature scalar TZ, the twist (as described by the normalized 
component Ug = UaO"", ), and the "plus" component of the outgoing shear 
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where 6°" and 0" are unit vectors in the 6 and directions. (The remaining components of 
the twist and shear vanish because of symmetry.) Because of axial and reflection symmetry, 
the range < ^ < 7r/2 suffices to display the full angular behavior. 

The first set of surface plots, Figs. 0-|n|, are for the mildly nonlinear case e = 10~^, 
corresponding to panel (c) in Fig. |^ in which the expansion wins the race. The evolution 
is traced from the starting time to the finish when the Bondi boundary forms. Figure |^ 
shows that the radius decreases fastest at the equator, in accord with the trouser-shaped 
horizon picture. Figure ^ shows that the expansion wins the race along a polar ray. Figure |^ 
shows that the curvature scalar TZ of the conformal horizon metric Hab increases from its 
unit sphere value 7?. = 2 in the region near the poles and decreases in the region near the 
equator. This is what would be expected for the conformal geometry of the Minkowski time 
slices t = const of of the collapsing prolate wave front which seeds the model. However, it is 
important to bear in mind that the quantities in Figs. |7-11 refer to the curved space affine 
time slices t = const. This is emphasized in Fig. ^ which shows the behavior of the twist, 
a quantity which would vanish for the Minkowski time slices of the fiat space wave front. 
The shear a-^-^uT, as depicted in Fig. |Tl|, is positive near the poles and negative near the 
equator, again as would be expected from a prolate Minkowskian wave front (according to 
our above polarization convention). Although the time dependence of all quantities in the 
model is analytic, there is an apparent "crease" at t ^= —80 in the surface plots of the 
Ricci scalar, twist and shear which results from a rapid change in the quantity A~^ = dt/dt 
governing the relative Minkowski and curved space affine times. For models with larger e, 
the angular dependences are qualitatively similar but the time dependence is more dramatic. 



VII. DISCUSSION 

Expressed now in terms of a black hole merger, this work has traced the horizon structure 
of a head-on black hole collision from the close approximation to the nonlinear regime. It 
has revealed dramatic time dependence in the intrinsic and extrinsic curvature properties 
of the horizon in the extreme nonlinear regime. The results suggest two classes of binary 
black hole space-times depending upon whether the crotch in the standard trouser picture 
is protected, in the sense that it lies inside a marginally anti-trapped surface on the horizon, 
or bare. Only in the bare case is it possible that the black holes are formed by either the 
collapse of matter or the implosion gravitational waves (see Fig. 0) originating in an initially 
nonsingular space-time. 

The results pave the way for an application of the PITT code to calculate the fully 
nonlinear wave forms emitted in the merger to ringdown phase. It remains to be seen in 
future work whether the dramatic time dependence in the merger stage is responsible for 
equally dramatic wave forms. 

While the numerical results presented in this paper are for the axially symmetric case, 
the codes are not restricted to any symmetry. It will be interesting to see how the results 
for a head-on collision are modified in the inspiral and merger of spinning black holes. 
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APPENDIX A: NUMERICAL INTEGRATION 

1. Integration of p 

As explained in Sec. |ITT|, the conformal model supplies the area coordinate r as well as 
the null data J, but for completeness we describe here the integration of the Raychaudhuri 
equation ( |2.5| ) that determines r from the null data in a more general setting. We integrate 
this equation in terms of the variable p = t/tm, where tm, defined in Eq. (|4.8| ), satisfies 
rM|A=o = and tm = 0. As a result, the evolution equation for p is identical to that for r, 
which in the spin-weighted form of Eq. ( p.27|) becomes 



p = -^{jj-k'). (Ai) 

The initial data consists of the values of p and p on the initial slice. Note that K = ^{JJ)/K. 
We put the equation in first-order form, 

p = H, il = -4Sp (A2) 

withS= {JJ-k^)/16. 

The advantage of the first-order form is that the pair of equations ( |A2| ) can be discretized 
to second-order accuracy using only two time levels, in the same footing as the other horizon 



evolution equations. The time integration stencil is the midpoint rule |23 



n+i _ n 1 

^ . ^ =-(H"+^ + H") (A3a) 

/\u 2 

-rrn+l _ Tin 

■ = -2^"+2(p"+l+p") (A3b) 



An 
which can be solved simultaneously for p""*"^ and H""*"^ to give 



p"(l - SAu^) + H"Am 
1 + SAu^ 



p"+i = ^^ I2 (A4a) 



H"- ^ mi-SAu^)^ASp^Au^ ^^^^^ 

1 + SAu 
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2. Integration of uj 



The time dependence of uj is determined by Eq. ( p.29|) , which we renormahze by factoring 
out r^ from both sides. On the left hand side, this is accomphshed by using the identity 

{r^u)=rlj{p'u), (A5) 

which holds on the horizon. Similarly, we re-express the first term on the right-hand side as 



r^a 






g 



rlp'^ 



(A6) 



Equation (|2.29|) then reduces to 
1 



(/^) 



4 



P2P + Pi ] p + Po, where 



(A7a) 



Po = -pQp (A7b) 

p^ = Qp + (jj - kK)Qp + {JK - KJ)Qp (A7c) 

P2 = JQJ - AKQK + 2kBj + 3 J5 J - 2Jm + 2Jdk + 2 JS J - 2/i 8 J - 27^37^. (A7d) 
We use the midpoint rule to integrate Eq. ( |A7a]) , i.e. the left-hand side is evaluated as 



(pM-^((p^ 



'u;r+' 



ip'^r) 



and p, J, etc. in Pq, Pi and P2 are evaluated at t"^^'^, e.g'. 



(A8) 



(A9) 



[Since uj does not enter in the right-hand side, this is a special case of the second order 
Runge-Kutta scheme ||2^, also used below in Eqs. ( [AllD and ( |A12D .] 



3. Integration of p\ 



The time dependence of r a, determined by Eq. (|2.3CI| ), is re-expressed in terms of p^\ 
using the ansatz r = tmP, which yields, for A = 0, 



{r'')^xu = ^M^PPM + (8MV,a - 2up)p,u - p'- 
Substitution into Eq. ( p. 301 ) then gives 



(AlO) 



%M^pp,Xu + (8MV,a - 2up)p,„ = ^ 



(g/f _ gj) ( j; + ^ j + A- (^8>^ + ^ - -^ 



P P 



Jtu^ + Kcju; - -n + p^ 



, (All) 



which we integrate using a second order Runge-Kutta scheme. 
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4. Integration of J\ 

After setting r = tmP, the evolution equation (|2.34D for Ja becomes 



SM'^p^J^Xu - up^j + (1 + K^) Scu = -8MV (pJa + P,\J) + 4MV ( JJ a + JJ,\ - 2 i^/s:,A) 

+ (1 + K^) (uj^ - 2 cj^ j - cj ( jBiT - K^J) - uj {KdJ - JdK) 

+J (jQlu -K{Qu; + Su) + Jlu^ - 2Kuju + 2 — [u^p + uQp) - 2 Juj-^\ (A12) 

which we again integrate using a second order Runge-Kutta scheme. 

APPENDIX B: THE CASE OF AXIAL SYMMETRY 



In axial symmetry the Sachs metric Eq. (2J.) can be written as 



ds^ = -{W - U^-ir~^)du^ - 2dud\ - 2Udu dx + r^ (— + 7 dipA . (Bl) 

Here W, U and r are functions of {u,\,x = cos 9). With the dyad choice associated with 
Eq. ( |3.14| ) the conformal horizon model gives 



7 = (i-x^)r, r=(^^J, j=-(i/r-r), (B2) 

with r a well-behaved function on the sphere. For the prolate spheroidal model considered 
here, F vanishes at the points t = r^ where the horizon pinches off and F = O^r"^). Thus, 
in the prolate case, it is useful for numerical purposes to use 7, as opposed to 7"^, as the 
variable to represent the metric on the two-sphere of constant u and A. The twist uj is related 
to the quantity U by 

uj = -^Vl^^U,x. (B3) 

The Einstein equations yield the following system of PDE's for propagating the metric 
variables along the horizon: 

Ou{r U^x) = 2 r n ^,x - 2 r r ^„ (B4) 

7 7 

a a 2 _ (^.x)^7 , ^,xx7 , {U,xf 1 U,Xxl . r^^ 7,^ ^,A 7,x , 7,xx f-n^. 

In the Schwarzschild case, r = tm = 2M — {Xu/4M), U = and 7 = (1 — x^). Thus, 
if the geometry is near Schwarzschild, all terms on the right-hand side are small except 
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7a;a; = —2. Ill Older to make Eq. (p5| ) numerically well behaved we subtract this term by 
introducing the auxiliary variable A = 9a (r^ — r|^). Then A satisfies 



a.A 



(rj^ ^ n^ ^ (Ml _ C^ ^ n^ _ fAvT. ^ /2j£ + 1 



The quantity r a is reconstructed as 



(B7) 



^A 



A + (r^ 



MJ,^ 



U 



2r 2r 

Special care is needed in order to write the right-hand sides of Eqs. ( P4|) and 



fB8) 



m 



manifestly regular form (for r 7^ 0) at the axis of symmetry {x = —1, 1) where 7 = (1 — a;^) F 
vanishes. We thus express 7,^/7 as T „/r and 7,x«/7 as 



7, 



2xT^ 



7 



r 



'l-x^)T 



(B9) 



Axial symmetry implies that r_„ vanishes at the poles so that the right-hand side of Eq. 
is regular, as can be seen by differentiating the conformal data in Eq. (p^ ) to obtain 



r,« 



/ 9f \ {r^ - re) 
[duj {i-ref 



(BIO) 



and noting that r^ = rg on the symmetry axis. Thus Eq. ( p4| ) is rendered explicitly regular 
on the symmetry axis. At the pinch-off points at r = 0, Eqs. (|B4|) - (|B6| ) are singular but 
the numerical performance near these points can be enhanced by noting that 7/r^ remains 
regular, by virtue of Eqs. (|3.11|) and ( p2D . 

The data at iS_ are initialized by the same prescription as for the general case discussed 



in Sec. [IV B| , which implies 



d^U. 



J_ 



2x 



:J- 



X^ 



(Bll) 



in accord with Eq. (|4.22|) 



A_ = n_ 1 



4MV 



in accord with Eq. ( |4.23| ) and 



dxl- 



0. 



(B12) 



(B13) 



in accord with Eq. (|4.23|) . 

We integrate Eqs. (P^), ( [BUD and ( [BT]) numerically to second order accuracy, using the 
same midpoint rule as for the general case in Appendix ^. 
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APPENDIX C: CONVERGENCE TESTS 

We have verified second order convergence of tlie liierarcliy of equations which provide 
horizon data for the 3D code. 

We have also checked second order convergence of the axially symmetric code, and we 
have used it to confirm the behavior, described in Sec. 0, of the sudden nonhnear plunge 
of the horizon radius to zero. Our aim is to use the axially symmetric code to generate an 
independent numerical solution against which to check the 3D code. 

Due to its lower computational requirements, the axisymmetric calculation is significantly 
more efficient than the 3D one. This will be particularly useful in the approach to the caustic, 
where the highest resolution is needed, and where the axisymmetric code will allow us to 
make more detailed studies of the behavior of the horizon. 
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FIGURES 




Source 



FIG. 1. The global binary black hole problem involves the event horizon TC~^ and future and 
past null infinity, here approximated by J'^ and J'~ . The sources which formed the black holes, 
e.g. collapsing matter or imploding gravitational waves, are indicated by the shaded regions. The 
sources are surrounded by a null world tube F, outside of which characteristic evolution supplies 
the exterior spacetime. 
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A = 



FIG. 2. The white hole horizon 7i is foliated by an affine parameter n, with S~ representing 
an early quasi-stationary slice. The affine parameter A along the outgoing null hypersurfaces Ju 
emanating from the foliation is chosen so that A = on H. The angular coordinates x^ are chosen 
to be constant along the light rays generating 7i and Ju- The null vectors /" and n" are used to 
project tensor fields into 7i. 
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FIG. 3. The ingoing null hypersurface normal to an spheroid Sq in Minkowski space pinches off 
in the future at points X, where two null rays cross, bounded by a caustic set C, where neighboring 
rays focus; and it expands in the past to a quasi-spherical slice S~. The foliation by Minkowski 
time t is is shown by horizontal slices. The conformal model induces a white hole afhne foliation 
u, indicated by dashed lines, whose upward bulge relative to the Minkowski foliation produces the 
fission into disjoint white holes depicted in the final n-slice shown. 
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FIG. 4. The "race" between the minimum values (over the sphere of rays) of the outward 
expansion and the radius of the horizon, for increasing eccentricities e = 10~^ (a), e = 10~^ (b), 
e = 10-5 (c), e = 10"'' (d), e = IQ-^ (e), e = lO'^ (f), all for mass M = 100 and to = -10. 
The expansion is plotted in terms of inf(0o;7T) with dashed lines and the radius in terms of 
inf (r/(2M)) with solid lines, both rescaled to unity at the starting time. For small e, the radius 
hovers near 2M longer than it takes the Bondi boundary to form and the fission remains hidden. 
For larger values of e, the fission occurs before the Bondi boundary, as the sequence shows. 
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FIG. 5. Embedding diagram for e = 10~^ to depict the shape of the surface, the scales being 
arbitrary. The darkened region indicates where the outward expansion has gone negative. The 
expansion first reaches zero at the poles. 
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FIG. 6. Embedding picture for e = 10~^. The radius reaches zero first at the equator, at 
the center of the picture where the individual white holes separate. The expansion is everywhere 
positive at that time. 
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FIG. 7. The radius decreases fastest at the equator, in accord with the trouser-shaped horizon 
picture. 
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FIG. 8. The minimum of the expansion occurs at a pole, hence the expansion wins the race 
along a polar ray. 
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FIG. 9. The curvature scalar TZ of the conformal horizon metric h^B increases from its unit 
sphere value 7?. = 2 in the region near the poles and decreases in the region near the equator. 
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FIG. 10. The behavior of the twist, a quantity which would vanish for the Minkowski time 
shces of the flat space wave front, shown here for e = 10~^. 
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FIG. 11. The shear gj^^qut ■, as depicted in Fig. 11, is positive near the poles and negative near 
the equator. 
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